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INTRODUCTION 


Even though the NASA Structural Analysis (NASTRAN) program 
is designed to solve numerous structural dynamic problems through 
the use of available rigid formats, an important class of problems, 
where the structures are spinning at a constant angular velocity, has 
been omitted. Rotating shafts, blades of spinning turbines, rotating 
linkages, and spin stabilized satellites are examples of problems 
falling within this class. These problems differ from the nonspinning 
structures in several significant ways. The accelerations of the masses 
in a nonrotating stationary frame are represented by the second deri- 
vatives with respect to time of the spatial variables. In the case of a 
structure spinning at a constant angular velocity, expressions for the 
accelerations of the discrete masses contain terms arising from the 
second derivatives of the spatial variables; in addition, they contain 
terms caused by coriolis accelerations, which are proportional to the 
velocities of the masses in the rotating frame. Finally, these expres- 
sions reflect the variations in steady-state centripetal accelerations 
caused by the small displacements of the masses in the rotating frame. 
The steady-state centrifugal forces set up the steady stresses that give 
rise to the geometric stiffness matrix. 

Since NASTRAN does not construct coriolis and centripetal 
acceleration matrices, and a centrifugal load vector due to spin about 
a selected point or about the mass center of the structure, a Fortran 
subroutine to construct these matrices is added in NASTRAN. The 
rigid translational degrees of freedom can be removed by using a 
transformation matrix T and its explicitly given inverse, T“V. These 
matrices are generated in the above Fortran subroutine and their 
explicit expressions are given in Appendix A. 
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The complex eigenvalue subroutine of NASTRAN does not 
measure up to the excellence it has shown in assembling the matrices. 
If the user desires, an option is available to write out the matrices 
generated by NASTRAN on a magnetic tape which, in turn, can be used 
as the input to another eigenvalue program. The probable advantages 
in using another eigenvalue program are that the user may be able to 
solve a larger problem within the available core and he may use a 
more efficient eigenvalue routine if one is available to him. If it is 
required, the user can write out certain information generated by 
NASTRAN on a magnetic tape unit using the subroutines OUTPUT2 
and WRTAPE used in this program. 


THEORETICAL DESCRIPTION 1 ' 2 


The equations of motion of a spinning structure are briefly 
derived here to show how they differ from those of a nonspinning struc- 
ture. The direct use of the Newton-Euler equations gives 


F s = m s A s 


T ^ 


*d H s 
dt 


( 1 ) 


for the s^h rigid body of a flexible appendage; where m s is the mass, 
A 8 is the absolute translational acceleration vector. F s and T s are 
the sum of the external and connection force and torque vectors, 
respectively. H s is the angular momentum vector and i denoted 
differentiation in the inertial frame of reference. 

For a rigid body of an appendage spinning nominally in the 
steady state with an angular velocity <o (fig. 1), the expression for 
acceleration is written as 
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where c is a vector representing the location of the mass center at 
time t with respect to its steady state position. u s and are vectors 
representing the displacement and small rotation, respectively, of the 
s th rigid body from its steady state configuration, r s is a vector 
representing the location of the rigid body of the appendage in its 
steady state configuration measured from the steady state mass center 
location. I_ s is the inertia dyadic of the s* b rigid body. Superscript b 
denotes differentiation in the reference frame b imbedded in the rigid 
body with the origin at the steady state mass center location. 

For zero spin (oj 

form 

A s 


T s = I s • — /3 s . (3-a) 

“ dt 2 ~ 


= 0), eq. (2) and (3) reduce to the familiar 
id 2 

= — — (R + c + r s + u s ) (2-a) 

dt 2 — — — 


In matrix notation the second term on the righ hand side of eq. (2), 
which is due to coriolis acceleration, gives rise to a skew- symmetric 
matrix; whereas, the third term, which is due to the centripetal accel- 
eration, yields a symmetric matrix. The fourth term in eq. (2) and (3) 
represents a steady state centripetal acceleration which describes the 
steady state configuration. Stretching forces, moments and rotations ob- 
tained therefrom, give rise to the second order geometric stiffness matrix. 
In the absence of angular acceleration, the fifth term of eq. (2) vanishes. 

If rotational dynamics are the primary concern, the effect of translation 
of the orbit is disregarded and the last term of eq. (2) also vanishes. 

The last two terms in eq. (3) will cancel each other if the inertia matrix 
is diagonal with all the terms having the same magnitude. In the com- 
puter program no such restriction is imposed on the I s matrix. 
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Conservation of linear momentum provides the relation 


c 


_L 

M 


n 

2 

s = l 


m s u s 


( 4 ) 


where M is the cumulative mass of all appendages and the central rigid 
body, and n is the total number of masses representing all of the 
appendages. 

Conservation of angular momentum is not imposed. As a result, 
the central rigid body is restricted against variations in rotations. 
Conservation of linear momentum permits the translation of the central 
rigid body, thus allowing the coupling of the vibrations of all the append- 
ages attached to the central rigid body. 


The set of equations representing the motion of all the rigid 
bodies in the appendage about the steady state configuration is obtained 
by substituting eq. (4) into eq. (2) and writing the resulting eq. (2) and 
(3), for all rigid bodies in matrix form; 

[M'j {u} + [G'] {u} + [K m + Ke + Kg] {u} = [F] . (5) 

The steady state equation in matrix form is: 

[Ke + K"'] {u} S = {P} S . (6) 

The use of eq. (4) eliminates the remaining translational rigid 
body degrees- of-freedom. As a result, the mass matrix M' is a 
symmetric non-diagonal matrix. Matrix G' is in general, a fully 
populated skew- symmetric matrix of coriolis acceleration terms. 
Matrix K"' is a fully populated non- symmetric matrix of centrifugal 
acceleration terms. Ke and Kg are elastic and geometric (differential) 
stiffness matrices, respectively, (and are obtained from the NASTRAN 
program) and {u} is the vector of generalized displacements about the 
steady state configuration. In the absence of spin, matrices G 1 , K"', 
and Kg will all be identically zero, and the eigenvalue problem reduces 
to the standard eigenvalue problem of a free-free structure or a canti- 
lever. { u} s is the vector of the steady state generalized displacements 
from the unstrained configuration {r} of the appendages. Since the 
steady state deformations {u} s are very small compared to the un- 
strained configuration {r}, it is assumed that the steady state config- 
uration is given by {r} instead of {r} + {u} s . The steady state force 
vector {p} s is used to obtain the geometric stiffness matrix. Kg . 
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Matrices M', G', and K" 1 have the following properties: 


[M'] 

= [M] 

[t] 

[G'l 

= [G] 

[t] 

[K"'] 

= [K"] 

[T] 


Relations (7) afford a transformation 

{y} = [t] {u> 

Substitution of transformation (8) into eq. (5) gives 

[M] {*y} + [G] {y} + [K" + [Ke + Kg] T" 1 ] {y} = 0 


(?) 

(8) 

(9) 


where 


u = Vector of displacements from the steady state config- 
uration of the nodal masses in spinning body frame. 

T = Transformation matrix relates the displacements of nodal 
masses in the body frame with the origin at steady state 
mass center to the displacements in another body frame 
obtained by translating the above frame to the instanta- 
neous mass center. In the absence of vibrations both 
above body frames coincide. If the axis of rotation and 
the origin of the body frame are both fixed in inertial 
space, T and T _1 become identity matrices, additionally. 

M = NASTRAN generated mass matrix 

G = Dummy module generated coriolis acceleration matrix 

K" = Dummy module generated centripetal acceleration 
matrix 

Pg = Dummy module generated steady state centrifugal 
force vector 

Kg = NASTRAN generated differential stiffness matrix using 
the above load vector Pg 

Ke = NASTRAN generated elastic stiffness matrix. 

Explicit forms of the above matrices are given in Appendix A. 
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DMAP DESCRIPTION 


The following information and options are made available through 
the input of vector WW with five elements using DMI* cards. The first 
of the two cards never changes for this program. WW(1), WW(2), 

WW(3) are the components of the spin vector in the body frame. Terms 
WW(4), WW(5) can take the values either 0. 0 or 1. 0. If WW(4) = 1.0 
the structure is spinning about the mass center of vehicle, and if 
WW(4) = 0. 0 the structure is spinning about a fixed point in the space. 
The calculation of matrices T and T , which removes the rigid body 
translational degrees of freedom, is performed if WW(5) = 1. 0. If 
WW( 5) = 0. 0, matrices T and T _1 are identity matrices which means 
that the structure is supported and does not have the rigid body trans- 
lational degrees of freedom. 

The following options can be exercised through the use of WW(4) 
and WW(5). 

Case I. WW(4) = WW(5) = 1. 0, GRID 1 constrained in all six 
directions. The structure is spinning about the vehicle mass center, 
and the rigid body translational degrees of freedom are removed. 

GRID No. 1 is connected by a rigid link to the mass center of the 
vehicle in the' steady state configuration and one or more appendages 
are cantilevered from GRID No. 1, GRID No. 1 should be constrained 
in all six directions by use of SPC cards or permanent SPC on GRID 
cards. 

Case II. WW(4) = WW(5) = 0. 0, GRID 1 constrained in all six 
directions. The structure is assumed to be spinning about a point 
(GRID 1) fixed in inertial space, e.g. , a spinning shaft with GRID 1 
at bearing. 

Case III. WW(4) = 1. 0 WW(5) = 0. 0, GRID 1 constrained in 
all six directions. Node No. 1 is rigidly connected to the steady state 
mass center which is fixed in inertial space. The structure is spinning 
about mass center. 

Case IV. WW(4) = 0. 0 WW(5) = 1.0, GRID 1 constrained in all 
six directions. The structure is spinning about GRID 1 with rigid body 
translational degrees of freedom removed. 


* Refer to NASTRAN User's Manual for definitions of card names 
used herein. 
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The DMAP sequence given in Appendix B solves eq. (9) and 
eigenvectors y^ (PHI in DMAP) thus obtained are transformed to u^ 
(PHID in DMAP) which are the eigenvectors of eq. (8). 

It is essential that two subcases are used in the case control 
deck as shown below for successful completion of the NASTRAN run. 

CASE CONTROL DECK 
TITLE 


SUBCASE 1 

DISPLACEMENT = ALL 


SUBCASE 2 

DSCOEFFICIENT = DEFAULT 
BEGIN BULK 


No provision for checkpoint is made since the time taken to 
assemble the matrices is just a fraction of the time taken to find a 
few eigenvalues. 


FUNCTIONAL MODULE PROGRAMING NOTES 


In writing a functional module for NASTRAN, the concept of 
open core should be employed even if the corresponding logic for an 
open core array is not used. This gives the generality and possibility 
of later expansions without having to alter the program extensively. 

This does not mean that the fixed dimensioned arrays cannot be used 
in NASTRAN functional modules. The details of the open core concept 
are given in Section 1, 5 of the Programer's Manual. Once the dimensions 
are set either by open core or by fixed locations, the next steps are 
either to retrieve the data (input blocks) to be used for further compu- 
tations or to store the computed data (output blocks) in a prescribed 
format within NASTRAN. The data as described in Section 2. 2 of the 
Programer's Manual may be in the form of a matrix, a table or bulk 
data cards. 

A matrix data are stored in two separate parts. One part 
constitutes the name of the matrix in alpha-numeric form (Header 
Information), The columns of the matrix are stored on random access 
peripheral equipments. 
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The second part is called the Trailer Information and is stored 
in FIAT which is an executive system table of NASTRAN. The first 
part is stored as a set of logical records; the first record is the 
Header information, and the second and subsequent, logical records 
until the end of file is reached are the columns of a matrix. The 
Trailer informations, which is the collection of the properties (size, 
real, complex, symmetric, etc.) of all the matrices used in NASTRAN 
are given in the Programer’s Manual but if new matrices are created, 
their Trailer information should be stored according to the instructions 
on page 2. 2-2 of the Programer's Manual. Either of the above two 
parts describing a matrix can be called, as shown later, without 
disturbing the other. 


Each of the matrices, whether constructed by NASTRAN or 
computed in a functional module and designated as an input or output 
block in a particular DMAP statement should be referred to by a file 
number. The numbering system of a file is standardized by NASTRAN 
as consisting of three digits. The first digit takes value 1 if it is an 
input block and 2 if .it is an output block; the second and third digits 
refer to matrix location in the string of input or output data blocks. 

For example, file number 102 in DMAP statement DUMMOD1 given in 
Appendix B refers to the second input block which is an unreduced mass 
matrix Mgg whereas 203 refers to the third output block which is the 
coriolis acceleration matrix G. 


In order to read the desired matrix, the following set of calls 
to subroutines listed below will unpack and read the matrix data. In 
each of the subroutines the file number for the appropriate matrix 
data block must be included in appropriate argument. 

CALL RDTRL 
CALL OPEN 
CALL FWDREC 
CALL UNPACK 


CALL UNPACK 
CALL CLOSE 

The subroutine RDTRL calls on the file number appearing in its argu- 
ment for the Trailer information. A call on RDTRL can also be made 
after calling OPEN, if desired. Subroutine OPEN opens the file to be 
read. FWDREC positions the requested file forward one logical record 
thereby skipping the first record in this particular example. If for 
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some reason two logical records need to be skipped, FWDREC is called 
twice. Each call to UNPACK allows the reading of one column (a logical 
record) of a matrix at a time. The call to UNPACK can be put within a 
DO-loop once the information on the number of columns of the matrix is 
obtained from Trailer information. After the reading of columns is 
completed the subroutine CLOSE is called to close the file as soon as 
practicable. 

If the data are in tabular form, instead of calling UNPACK, call 
READ to read the data. Each call to READ reads one logical record 
of the data. The programer's manual should be consulted for structure 
of the record read. The call to READ can be put either within a DO-loop 
once the information on the number of records is obtained from Trailer 
information or within an unending DO-loop in which case, when the end 
of file is reached, the transfer will be made to a statement number 
appearing in the argument of READ. The set of calls is shown as 

CALL RDTRL 
CALL OPEN 
CALL FWDREC 
CALL READ 


CALL READ 
CALL CLOSE 

If the data on Bulk Data Cards are desired to be retrieved the 
following set of calls to subroutines should be employed. 

CALL PRELOC 
CALL LOCATE 
CALL READ 


CALL READ 
CALL CLOSE 

Subroutine PRELOC locates the file on which the bulk data card 
images are stored and LOCATE locates the desired number of cards 
in the file. 

In DMAP statement DUMMOD1 file number 101 which is GEOM1, 
.contains the geometric information from Bulk Data Cards. This file 
number is called in PRELOC subroutine. For each bulk data card 
to be read, subroutine READ should be called. 
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To pack the matrices calculated in a subroutine and appearing 
as output data block the following set of calls to subroutines is required. 

CALL OPEN 

CALL FNAME (finds Header information) 

Prepare Trailer information (e. g. M(l), M(2),--M(7) ) 
according to instructions appearing in Section 2.2-2 
of the Programed s Manual. 

CALL WRITE (writes Header information) 

Perform computations. 

CALL PACK (----, ----, WRITE, M) 


CALL PACK 

CALL WRTTRL (M(l)) 

CALL CLOSE 

Subroutine OPEN opens the file to be written on and FNAME 
finds and stores the Header information as appearing in DMAP sub- 
routine (e. g. , if file number 204 is referred in OPEN, FNAME will 
go to the fourth output block of DUMMOD1 which is matrix AA and 
store AA as the Header information). Next prepare the Trailer infor- 
mation in the vector M according to the instructions in Programer's 
Manual with the following exceptions. • On page 2. 2-2 of the Programer's 
Manual M has the dimension 6 which is an error, it should be dimen- 
sioned M(7) and M(l) = File Number 
M(2) = 0 

M(3) = M(2) of the Programer's Manual 


M(7) = M(6) of the Programer's Manual. 

Note that M(2) is set equal to zero and not one as implied in the 
Programer's Manual. If M(2) - 1, the information on the number of 
columns stored in Trailer information will show one more than the 
value desired. Hence when this Trailer information is used to read 
the columns, the READ will try to take it past the end of file resulting 
in fatal error. 

Subroutine WRITE writes the Header information and then one 
column at a time is packed by subroutine PACK. Call WRTRL to write 
Trailer information and then call subroutine CLOSE. 
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If it is desired when developing functional modules, the number 
of input, output blocks and parameters can be altered by altering the 
information in MPL (subroutine XMPLBD). In the development of this 
program, four functional modules are written. In addition several 
other functional modules for other NASA projects at MSFC have also 
been written. All these functional modules are given dummy names as 
given in User's Manual. Because the procedure given in NASjTRAN's 
program manual are either incomplete or in error, attempts to add 
new functional modules with unique name and unique MPL (Module 
Properties List) have not been met any degree of success. 

It is acknowledged that most of the information presented in this 
section may be found throughout the Programer's Manual. However, 
it will take a considerable time to assemble and use. Here we have 
presented the information we collected through trial and error and 
months of diligent work by two expert programers. We present it 
with the hopes that someone wishing to write their new functional 
modules will not have to encounter the same difficulties. 
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APPENDIX A 
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The inverse of this matrix is expressed in the following form, rather than 
inverting [T] by some matrix inversion technique. 


r 

\ 


10 


0 1 + v ’* 0 

20 , 


0 14 ^- 

S 30 


o <r o 

ooo 
o o o 



o o o 
ooo 
0 0 0 



0 0 
0 0 
0 0 


IT] 


1 



( ooo 
0 0 0 
ooo 


10 0 ooo ooo 


0 0 0 0 0 


0 10 ooo 


000 ... 000 00 


0 0 1 ooo ooo 
ooo 

ooo 

0 0 0 t 




0 

0 



0 

0 

0 


0 



0 0 
0 0 

0 0 

0 0 


I 

ooo 

t 


0 0 0 0 0 


ooo 


0 0 0 0 0 


ooo 


0 0 0 0 0 



ooo 

ooo 

ooo 

ooo 

ooo 

ooo 


b 0 

10 

0 l«hr~ 
20 


30 

0 

0 

0 



6nx 6n 


213 



where 
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mi 

n 

rM 0 + 2 mJ 

j = l 
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if the motion of mass at node i in kth 
direction is not constrained to be 
zero 


if the motion of mass at node i in kth 
direction is constrained to be zero 


S ko = S / 

k 


where m 1 is the mass at the ith node point of the total of 'n' nodes and 
ijk is the moment of inertia of the rigid body at the ith node, M 0 is the 
mass of the central rigid body and {&>} is the spin vector. 

Matrices [G] , [K n ] , [T] , and [T -1 ] are the square matrices of 
the dimension 6n x 6n. Rows and columns corresponding to the degree- 
of-freedom which are either constrained to be zero or have no mass 
should be removed. This will reduce the above matrices to NxN where 
N is the total degree-of-freedom of the problem. 
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APPENDIX B 
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123/V,N,.N0SIMP/C,NiO/V<N,N0GENL/V,N,GENEL S 
N0GENL,IN0SIMP,GENEL $ 

ERR0R1/.N0SIMP $ 

0GPST/GENEL * 

CSTMjHPT, ECPT ,GPCT 1 01 T /KGCX ,.»GPST/V»N»N0GENL/V,N» N0K4GG * 
CS:TM,MPT,ECPT,GPCTiOIT7MGGi/V,Y,WTMASS=l.O/V,N,N0MGG/V,N,N08G3/ 
V,Y,C0UPMASS=~1 $ 

N0MGG 4 
LBL1.GR0PNT $ 

ERR0R4/N0MGG $ 

BGPOT,CSTMIEQEXIN,HGG/0GPWG/V,Y,GRDPNT=-L/V,Y,WT«ASS * 

0GPWG, /,.» , //V»N»CARDN0 * 

CARDN0 $ 

LBL1 $ 

KGGX , KGG/N0GENL $ 

LBL11 , N0GENL S 

GE I , KGGX/KGG/V, N ,LUSET/V ,NtN0GENL/V,N»N0SIMP $ 

LBLll * 

//CiN,MPY/V,N,NSKIP/C,N,0/C,.N,0 i 

C ASECC, GE0M4, EQEXIN, SIL,GPDT/RG,YS^USET/V,N, LUSET 7V,N,MPCF1/V, 
N,MPCF2/V,N,SINGLE/V,N<0MIT/V,N,REACT/V,N,NSKIP/V,N,REPEAT/V, 
N,N0SET/V*N*N0L/V,N,.N0A $ 

MPCFl,MPCF2,SINGLEk0MIT„REACT,NSKIP,REPEAT,N0SET, N0L.N0A $ 
ERR0RS*N0L * 

GM/MPCF1/G0, K00B ,L00 , U00 ,P0 , U00V, RUB V/0M I T/PS »KFS »KSS ,QG/ 
SINGLE/ UB 00 V/0MIT/YBS, P BS, KBFS KBSS,KDFS« KDSS/SI M3LE $ 
KGG,KNN/MPCF1/MGG,MMN/MPCF1 $ 

LBL4D, REACT $ 

ERR0R2 $ 

L8L4D % 

LBL4,GENEL * 

GPL,GP3T,USET,SIL/0GPST $ 

0GPST t i ft t // V »N'» CARON0 $ 

CARDN0 $ 

L8L4 * 

LBL2.MPCF2 $ 

USET , RG/GM > 

USET »GH».K:GG»MGG» ,/KNN,MNN, I $ 

LBL2 t 

KNN,KFF/SINGLE/MNN*MFF/SINGLE $ 

LBL3, SINGLE $ 
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53 sc ex 

54 LABEL 

55 EQUIV 

56 C0NO 

57 SMP1 

58 LABEL 

59 DUMM0D1 

60 EQUIV 

61 R.B.MG2 

62 C0NO 

63 SSG2 

64 LABEL 

65 SSG3 

66 C0ND 

67 MATGPR 

68 MATGPR 

69 LABEL 

70 SDR1 

71 DUMM0D2 

72 SAVE 

73 MATPRN 

74 C0NO 

75 DSMG1 

76 SAVE 

77 ADD 

78 EQUIV 

79 C0ND 

80 MCE 2 

81 LABEL 

82 EQUIV 

83 C 0 NO 

84 SC El 

85 LABEL 

86 EQUIV 

87 C0ND 

88 SMPt 

89 LABEL 

90 LABEL 

91 EQUIV 

92 MATPRN 

93 MPYAD 

94 MPYAD 

95 ADD 

96 DPD 


97 0UTPUT2 

98 0UTPUT2 

99 CEAD 

100 SAVE 

101 0FP 

102 SAVE 

103 C0ND 

104 MPYAD 

105 MATPRN 

106 M0DG 

107 VCR 

108 SAVE 


USET 7 KNNVMN^ ** /K FFi KFSfK $$ i MFFV * $ 

LBL3 $ 

KFF9 K AA/0M IT /MFF »MAA/0M I T '$ 

LBL5 1 0.MIT $ 

USET » KFFf MFF f f /G0, KAA , K0 08 i L00^ U00 , MAA, M008 ,M0AB, , $ 

LBL5 $ 

GE0M1,MGG,BGPDT, WV^USET, , , /KP,PG ,G:» AA,T, TI ,RPP,MD0F/ $ 

PG i PL/N0S ET $ 

KAA/LLL tULL $ 

LBLIO ,N05ET $ 

USET t GM >Y$ J KFS LG0t TPG/ fP0, PS , PL $ 

LBLIO S 

LLL,ULL ? ,KiAA,PLa00^U00fK008»P0/ULV,U00V,RULV,RU0V/V,N,0MIT/ 

V, V , I RE S=-l $ 

LBL9, IRES $ 

GPL tUSETrSIL kRUL V/ /C > NV L $ 

GPL » USE T I L ,RU0V/ /C » N $0 $ 

LBL9 $ 

USET t PG *AJLV»U00V fY5**G0 iGMfPS ,KFS »K5S,/USV ,PGG ,QG/ C , N, 1/C,N,DS0$ 
UGV , , , i , / UGVX t m ? m t/VfYi I UGV=0 $ 

IUGV % 

UGVX » Wv// $ 

LBBtIUGV $ 

CA3ECC »GPTT » 4 SIL, EOT ,UGV, C STM , MPT ,ECPT »GPCT,DI T/KPSG/ 

V , N* DSC0SET $ 

DSCBS6T £ 

KPGG , KGGX/KDGG/C ,Y $ALPHA=M U 0, 0.0 ) /C , Y, 8ETA= M . 0, 0. 0 ) $ 

KDGG, KDNN/MPCF2 $ 

LBL2D,MPCF2 $ 

USET,GM*KDGG,, */KDNH, ,* $ 

LBL2D $ 

KDNN, KDFF/S ENGLE $ 

LBL3D, SINGLE $ 

USET » KDNN, \ ,/KDF F, KDFS IKDSS , , , $ 

LBL3D % 

KDFF, KDAA/0MIT $ 

LBL5D, 0MIT $ 

USET » KDFF , I ,7000 ,KDA A ,K00B1 *L001 , U00 1 , , , , , $ 

LBL5D $ 

LBB $ 

KAAi KDAA/ IUGV/G0 ,G00/ I UGV $ 

KAA,G0V *, // $ 

KDAA »TI k/KSUM/C , N» 0/C ,N, l/CrN,0/C »N, 1 $ 
KAA,TU/KSUM2/C,N,0/C,N,1/C,N,0/C,N,1 $ 

KSUM, KP/K SUM 1/C, Y, ALPHAS { l*0f0*0)/C,Y,BETA=(l. 0,3*0) $ 
DYNAMICS^GPL ,SIL ,USET/GPLD t SI LD, USETD ,TFP00L , , , , , » EED,EQDYN/ 

V, N, LUSET /V, N ,LUS£TD/V IN , N0TFL/ V , N , N0DLT/V, N , N0PS DL/ 
V f NiN0FRL/VrN,N0NLFT/VIN,N0TRL/V,-N,N0EED/C,N,123/V f N,N0UE $ 

RPP, ND0F , T-I , +/ /C ,Y VPl-^l/CLY , P2 = -l $ 

MAA, KSUM2 , KSUMVKP, G/ /C I Y » Pl*-1/C , Y, P2~-l $ 

KSUM1 , G ,.M AA, EED, CASECC/PHI $CLAMA»0CEIGS/V,Nt EIGVS $ 

EIGVS $ 

0CE I GS V CL AM A *. , f , // V *N , C ARDN0 .$ 

CARDN0 $ 

LBL 16, EIGVS $ 

II /PHI //PHID/C IN ,0/C , NI 1/Cf N, Q/C »N, 1 * 

PHID, PHI, , < // $ 

AA^PHID/ /Cl N*— 1 $ 

CASECC,'EQDYM^USETD;PHIOVCLAMA,^/0PHID,/C,N,CEIGN/C,N,DIRECT/ 

C * N* 0/V ,.N ,N0D/V, N, N0P/C , N,0 $ 

N0D* N0P $ 
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109 C0ND 
no 0 fp 

111 SAVE 

112 LABEL 

113 C0ND 
114EQUIV 

115 C0NO 

116 SDRl 

117 LABEL 

118 SDR2 

119 0FP 

120 SAVE 

121 LABEL 

122 JUMP 

123 LABEL 

124 PRTPARM 

125 LA8EL 

126 PRTPARM 

127 LABEL 

128 PRTPARM 

129 LABEL 

130 PRTPARM 

131 LABEL 

132 END 


LBL15.N0D t, 

0PHID,,k»»//V,N,CARON0 $ 

CARDN0 $ 

LBL15 $ 

LBL16»N0P i 
PH1D*CPHIP/N0A $ 

LBL17 rN0A $ 

USETO, ^PHlDj , ,G00,GM, ,KFS, V/CPHIP,yQPC/C» N,i/C i N, DYNAMICS $ 
L8L17 i 

CASBCC^CST.M»MPT»OIT»-EQDYN»SILDi»i » CL AM A » QPC » CPHI P » EST t / » 0QPC1 » 
0CPHI P jl0ESC 1 »0EFC1 4 / C »N *CE IG $ 

0CPHIP/0QPC1,0EFCH0ESC1, , // V ,N, CARDN0 $ 

CARDN0 $ 

LBL16 * 

FINIS $ 

ERR0R1 $ 

//C * N »-l/C < N >.M0OES$ 

CDDftD O * 

/ /C,N,-2/C IN *D IF FS T IF $ 

ERR0R4 S 

//C t N , -^/CiN *D I FFS7 IF $ 

ERR0RS $ 

//€* N»-5/C|N ►OIFFSTLF $ 

FINIS $ 

* 
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Description of DMAP Operations for Eigenvalue Analysis of Rotating Structures 

2. GP1 generates coordinate system transformation matrices, table of 
grid point locations, and tables for relating internal and external grid 
point numbers. 

4. GP2 generates Element Connection Table with internal indices. 

5. PLTSET transforms user input into a form used to drive structure 
plotter. 

7, PRTMSG prints error messages associated with structure plotter. 

10. Go to DMAP No. 14 if no undeformed structure plot request. 

11. PLOT generates all requested undeformed structure plots. 

13. PRTMSG prints plotter data and engineering data for each undeformed 
plot generated. 

15. GP3 generates Grid Point Temperature Table. 

16. TA1 generates element tables for use in matrix assembly and stress 
recovery. 

18. Go to DMAP No. 123 and print error message if there are no structural 
elements. 

20. SMA1 generates stiffness matrix [ K X ] and Grid Point Singularity 

Table. * gg 

21. SMA2 generates mass matrix [ M ] . 

23. Go to DMAP No. 28 if no weight and balance request. 

24. Go to DMAP No. 127 and print error message if no mass matrix exists. 

2 5. GPWG generates weight and balance information. 

26. OFP formats weight and balance information and places it on the system 
output file for printing. 
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29. Equivalence [ K ] to [ K ] if no general elements. 

gg gg 

30. Go to DMAP No. 32 if no general elements. 

X 

31. SMA3 adds general elements to f K ] to obtain stiffness matrix f K 1 

L gg gg 

34. GP4 generates flags defining members of various displacement sets 
(USET), forms multipoint constraint equations [ R ] ju | = 0 and 
forms enforced displacement vector |Y g J, g g ^ 

36. Go to DMAP No, 129 and print error message if no independent 
degrees of freedom are defined. 

38. Equivalence f K 1 to f K 1 and f M ] to f M ] if no multipoint 
constraints/ gg ' ' gg 


39. Go to DMAP No. 41 if no free- body supports supplied. 

40. Go to DMAP No. 125 and print error message if free- body supports 
are present. 

42. Go to DMAP No. 46 if general elements present. 

43. GPSP determines if possible grid point singularities remain. 

44. OFP formats table of possible grid point singularities and places it 
on the system output file for printing. 

47. Go to DMAP No. 50 if MCE1 and MCE2 have already been executed for 
current set of multipoint constraints. 

48. MCE1 partitions multipoint constraint equations f R ] = [ R I R ] 

and solves for multipoint constraint transformation matrix m| 

[ G ] = -[ R ] “1[ R ]. 

1 m L m J L n J 

49. MCE2 partitions stiffness and mass matrices 


[ K ] = 

1 gg 


*» 

| I 



i 

K 

nn 

' K 

nm 

+ “ 
1 K 

and 

u 

S' 5 " 

M M 

nn nm 

_ 

K 


M 

M 

mn 

mm 



mn 

mm 
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and performs matrix reductions 

[ K ] = [ K ] + [ G T ] [ K ] + f K T ] f G ] 
L nn J L nn J 1 rrr L mn J • mn J 1 m 


f f G 1 f K 1 r G ] and 
L m J L 1 m J 


r M 1 = f Ml + [ G T ] [ M ] + [ M T ] [ G ] 

L i nn J L m J L mn J L mn L m J 


+ [ G ] [ M ] [ G ] . 

1 m L mm L m 

51. Equivalence f K ] to f K .] and f M ] to f Ml if no single- 

1 nn L ff L nn 1 ff 

point constraints. 

52. Go to DMAP No. 54 if no single- point constraints. 

53. SCE1 partitions out single-point constraints. 


K.. i K 
ff fs 


M ff l M £ S 

f M nJ ~ M 

sf I ss 


[ K nJ= K~ + K- 
sf I ss 


55. Equivalence [ K ] to [ K ] and [ M ] to [ M ] if no omitted 

xx aa lx aa 

coordinates. 

56. Go to DMAP No. 58 if no omitted coordinates. 

57. SMP1 partitions constrained stiffness and mass matrices 



K 

1 K 



I M 

[ K ff ] = 

aa 


and 

* 

H-h 

i i 

!! 

aa ao 

— — 

K 

l K 


M { M 


oa 

i OO 



oa | oo 

• i — 


solves for transformation matrix [ G 1 = -f K 1 TK 1, 

L o J 1 oo 1 oa J 

and performs matrix reductions [K ] = [K ]+[ ] [ G ] 

aa aa oa o 

and [ M aa ] = [ M aa ] + [ [ G ]' + [ [ M^J 


+ [ G ] [ M ] r G ] . 

L o J L oo J 1 o 
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59. Dummy module DUMMOD1 constructs coriolis acceleration matrix 
[ G] , centripetal acceleration matrix [ K ], transformation 
matrix [ T] and its inverse [ T] \ and centrifugal load vector 
{P J . Rows and columns corresponding to the degrees of freedom 
constrained to be zero or have no mass have been removed from 
[ G] , [ K™ ] , f T] , and [ T] Centrifugal load vector j P l 

is in g- set and is reduced in the following D- MAP statements, ^ 


60. Equivalence jP^j to jP^j if no constraints applied, 

6 1. RMBG2 decomposes constrained stiffness matrix [ K = f L-^] [ U^] 

62. Go to DMAP No. 64 if no constraints applied. 


63. SSG2 applies constraints to static load vectors 



jp l = 

P 

• + 

[ G ] 

\ p i , 

1 nj 

l n 


L m 

1 ml 



o 


l p £ } = j 

ft} 

> [ 

K £ 8 ] | Y s' 

1 ’ 

l p i) ■ i 

! P a) 

+ [ 

G oJ | P o| 

► o 


65. SSG3 solves for displacements of independent coordinates 




t 


solves for displacements of omitted coordinates 


1-°} 

= [K ] - 1 

(p }, 

l o 1 

L oo 

1 or 


calculates residual vector (RULV) and residual vector error ratio 
for independent coordinates 


l 6 P il = l p i) 



6 1 


■ h T l l 5P il 
l p i T ) W 
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and calculates residual vector (RUOY) and residual vector error 
ratio for omitted coordinates 


I 5 p 1 


{ p 1 


K 

{u°l 

l o) 


l o' 


oo 

l oJ 


f 


l"ol 

ii 

l 6P ol 

K: 


L°1 

l U ol 



66. Go to DMAP No. 69 if residual vector is not to be printed. 

67. Print residual vector for independent coordinates (RULV). 

68. Print residual vector for omitted coordinates (RUOV), 

70. SDR1 recovers dependent displacements 


l U o) 

= 

|°o] | U ll + 1 

O 0 


u 

a 

u 

o 

s 1 

"fi ■ |! 

.B ■ ^ 
II 

| U n) , 

{V} 

= 

[ G m] P„l • 

u 

n _ 
u 

mi 

l"gl 


and recovers single-point forces of constraint 



71. DUMMOD2 checks if vector {Ugj is a null vector. IUGV = - 1 if {ug} 
is null (geometric stiffness matrix KDGG is also a null matrix) 
otherwise IUGY = 0. 

74. Go to DMAP No. 90 if IUGV = -1, 

75. DSMG1 generates differential stiffness matrix [KggP] , 

ADD elastic and geometric stiffness matrices in g-set 

[ K gg] + f K ggl ^ t K gg^ ' 


77. 



78. Equivalence [Kgg] to [K d n ] if no multipoint constraints. 

79. Go to DMAP No. 81 if no multipoint constraints. 


80. MCE2 partitions differential stiffness matrix 

K d | K d 

[K* ] = 

. I d 


mn ' mm 


d. " ci T ci 

and performs matrix reduction [K ] = [K ] + TC 1 [K 1 
r nn L nn J L L mn J 


+ [K d ] [G ] + [G T ] [K d ] [G ] . 
mn m m mm L m 


d d 

82. Equivalence [K ] to [K„,.] if no single-point constraints. 

nn ff 

83. Go to DMAP No. 85 if no single-point constraints. 

84. SCE1 partitions out single-point constraints. 


[K U ] = + 

nn , . , 

K d | K d 
sf I ss ^ 

d cl 

86. Equivalence [K f , ] to [K ] if no omitted coordinates. 

ii aa 

87. Go to DMAP No. 89 if no omitted coordinates. 

88. SMP1 partitions constrained stiffness matrix 


= 


"d | d 
K K 

aa | ao 


K | K | 
oa I oo | 


did 

solves for transformation matrix [G 00 ] = - [K ] [K ] , 
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d — d T d 

and performs matrix reductions [K ] = [K 1 + [K 1 [G 1 

aa aa oa oo 

d 

91. Equivalence [K ] to [K ] and [G ] to [G ] if geometric 

3-ct 3.E O OO 

stiffness matrix [K^ ] is a null matrix. 

gg 

d 1 

93. Multiplies the matrices [K ] [T] = [KSUM ] . 

3.3. 


94. Multiplies the matrices [K ] [T] ■= [KSUM2 ] . 

33 ' 


95. Adds matrix KSUM and the centripetal acceleration matrix [K m ] . 

96. DPD generates flags defining members of various displacement sets 
used in dynamic analysis (USETD), tables relating internal and 
external grid point numbers, including extra points introduced for 
dynamic analysis, and prepared Transfer Function Pool and Eigen- 
value Extraction Data. 

97. Matrices [Rp ] , [NDOF] and [T ]" * are output on magnetic tape. 

[Rp] is (n x 4) matrix where n = no. of grid points. First three 
columns represent the coordinates of grid points in basic coordinates 
and fourth column stores the mass. data at grid points. 

[NDOF] is (3 x n) matrix. Value of 1. 5 is written if the translational 
D.O. F. at a grid point is not constrained by SPC, MPC, OMIT or 
permanent SPC on GRID cards. Otherwise it is written 0. 0. 

98. Matrices [M aa ] , [K aa ] , [K aa J, [K'"] , [G ] are output on 
magnetic tape, 

[K aa ] is the reduced elastic stiffness matrix 

[K^ a ] reduced (elastic + geometric) stiffness matrix 

[K ,n ] is the reduced centripetal acceleration matrix 

[G ] is the reduced coriolis acceleration matrix. 

99. CEAD extracts complex eigenvalues from the equation 

C M ddP 2 + B ddP + K dd] l u dl = 0 
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101. 

103 . 

104 . 


105 . 

106 ,. 

107 . 

109 . 

110 . 

113 . 

114 . 

115 . 


and normalizes eigenvectors according to one of the following user 
requests: 

(1) Unit magnitude of selected coordinate 

(2) Unit magnitude of largest component. 

OFP formats the summary of complex eigenvalues and summary of 
eigenvalue extraction information and places them on the system 
output file for printing. 

Go to DMAP No. 121 if no eigenvalues found. 

{0}, the eigenvector of 

[Mp 2 + Gp + [K ,n + [K e + K g ] T' 1 ]] {0} = 0 

is given by complex eigenvalue analysis step #91. 

{0d} » the eigenvector of 

[MTp 2 + GTp + K ,n T + K e + K g ]{0 d } = 0 

is obtained in this step {0 d } = [ T] “ ^ {0} . 

Eigenvectors {0 d } and {0} are printed. 

[ 5 ^ a (3 x 3) matrix for each of the eigenvector { 0 d } is 
constructed and printed. 

YDR prepares eigenvectors for output, using only the independent 
degrees of freedom. 

Go to DMAP No. 112 if no output request for the independent 
degrees of freedom. 

OFP formats the eigenvectors for independent degrees of freedom 
and places them on the system output file for printing. 

Go to DMAP No. 121 is no output request involving dependent 
degrees of freedom or forces and stresses. 

Equivalence {0 d }to{0p} if no constraints applied. 

Go to DMAP No. 117 if no constraints applied. 
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116. SDR1 recovers dependent components of eigenvectors 


<Ad 


= [ G oo H0 d }> j ~cf o '\ = {<P f + ^ e} ’ 


j^>f+ <£ e 

r-v— 


i - {.rf* 4" } = [G ] {0 + 0 } > 

f LT n ^e J L r m J L L n e J 


0 +0 

r e 


0 


m 


| ■ < v 


and recovers single-point forces of constraint 

{q g } ~ [ K fg] {<j>f } * 

118. SDR2 calculates element forces and stresses (OESC1, OEFC1) 
and prepares eigenvectors and single-point forces of constraint 
for output (OCPHIP, OQPC1). 

119. OFP formats tables prepared by SDR2 and places them on the 
system output file for printing. 

122. Go to DMAP No. 131 and make normal exit. 

124. Normal mode analysis error message No. 1 - Mass matrix 
required for real eigenvalue analysis. 

126. Static analysis with differential stiffness error message No. 2 - 
Free body support not allowed. 

128. Static analysis with differential stiffness error message No. 4 - 
Mass matrix required for weight and balance calculations. 

130. Static analysis with differential stiffness error message No. 5 - 
No independent degrees of freedom have been defined. 
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nonooonoonono oao 


SUBR0UTINE DUM0O1 

INTEGER TYPING TYP0UT % TOGF » SYSBUF , C0RS2 
REAL M, MU { 300 ) 

EXTERNAL WRITE 
EXTERNAL READ 

DIMENSION HEADt'SH HEAD2(5), HEAD3(5), HEA>4<5), ISUD(2) 

D 4 IC 0NM ( 2} , A ( 8 ) , B(13), ft(4) t . RG(4}* W{5), M(3G0) 

D TOGF ( 300 ) , NDEGF ( 6 9 300 ) « NFRE(300), XKPU) 


0 # AA (3 1 3.00') » 

0 ^RHEAD ( 5) 

DIMENSION XI 11 1 300 ) 


RP { 300 *3.) t* MM f 7 ) 


0 

0 


C0MM0N 

C0MM0N 

C0MM0N 


* XI 211 300) , 
^XMGGC 300) , 
/DUM0PN/ XKP 
/PACK X/ TYPIN, 


X 1.22 C 300) ,♦ 
O0F(6 f 300) 


XI 31(300)* XI 32 ( 300) » XI33C3D1) 


TYP0UT, 


/UNPAKX/ I TYPE fr JJ( N, 


II, Ml, 
JINCR 


I NCR 


DATA N A M 1 , NAM2, NAM3 V NAM4 / 201'* 202, 203, 204 / 


DM100100 

DM100200 

DM100330 

DM100400 

DM100500 

DM100600 

DM100700 

DM100800 

DM100900 

DM1010D0 

DM101130 

DM1012D0 

DM101300 

DM101400 

DM101500 

DM101600 


C0MM0N /SYSTEM/ SYSBUF* 

0UTAPE 

DM101700 

C0MM0N XXX 



DM101800 

EQUIVALENCE (Alan ID) 



DM101900 

1 ,IAD7)% I A.) 



DM102000 

2 ,(&12U IG) 



DM102 IOO 
DM102200 

0PEN C0RE ARRAY 



DM102300 

OM102400 

EQUIVALENCE (XKPt2000), 

NDEGF) 

DM102500 

E , (XKP (3800 ) , 

M 

} 

DM102600 

E » l XKP ( 41 00 ) » 

NFRE 

) 

OM102700 

E , ( XKP ( 4400 ) , 

TDGF 

) 

DM1G2800 

E , (XKP (4700), 

A 

) 

DM102900 

E , (XKP (4708), 

B 

) 

DM103000 

E ,(XKPt4721>* 

R 

1 

DM103100 

E .* (XKP ( 4725 )> 

HEAD 

) 

0M103200 

E , ( XKP ( 47 30 ) * 

HEAD1) 

DM103300 

E , ( XKP { 4735 ) , 

HEAD2) 

DM103400 

E , < XKP (4740), 

HEAD3) 

DMI 03500 

E *( XKP (4745), 

HEAD*) 

DM103630 

E * ( XKP ( 4750), 

AA 

) 

DM1037Q0 

EQUIVALENCE (XKP(5650), 

RP 

) 

DM103800 

E ,( XKP (6550), 

MU 

) 

DM1Q3900 

E * ( XKP { 68 50 ) , 

XIII 


DM104030 

E , { XKP (7150), 

XI 21 

) 

DM104100 

E , ( XKP ( 7450 ) , 

XI 22 

) 

DMI 04200 

E , ( XKP (7750 ) , 

X 1 3 1 

) 

DM 104300 

E , ( XKP (80 50), 

XI 32 

) 

DM104400 

€ , (XKP ( 8350) * 

XI33 

> 

DM104500 

E , (NDEGF* O0F ) 


DM104600 


DM1C4730 


ir*0 *•**»•#****,** #»*****«*r*****Ht**»**iK«^»»*» ************* ********** 1048 30 


IFL1101) REFERS T0 GE0MI-THE FIRST INPUT DATA BL0CK 0F DMAP 
STATEMENT DUMM0D1* 

IC0NM REFERS T0 MGG-THE SEC0ND INPUT DATA BL0CK 3F DUMM0D1. 
IBG REFERS T0 BGPOT-THE THIRD INPUT DATA BL0CK 0F OUMM0D1 • 

I WW ( 104 ) REFERS T0 WW-THE BL0CK C0NTAINING W-$ AND R-T FLAGS* 
NAM5 REFERS T0 USET-THE FIFTH INPUT DATA BL0CK 0F DUMM0D1. 

THE GENERAL PR0CEDURE F0R READING A DATA BL0CK IS-FIRST, 

0PEN THE FILE CONTAINING THE BL0CIC SUCH AS IC0NM,IB3 ETC . 

THEN SKIP THE HEADER RECORD BY CALLING FWDREC • THE NEXT STEP 
IS T0 DETERMINE IF THE DATA IS A MATRIX EG. MGG • F0R MATRICES, 


DM 104900 
DM105000 
DM105100 
DM105200 
DM105300 
DM 105400 
DM105500 
DM105600 
DM105700 
DM105800 
DM105900 
DM106000 
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THE NEXT STEP IS T0 CALL UNPACK. EACH CALL T0 UNPAC< BRINGS INT0 
MEM0RY 0NE COLUMN 0E THE MATRIX. F0R N0N-MATR I X INPUT* THE RECORD 
STRUCTURE MUST BE L00KED UP IN THE NASTRAN PROGRAMME R-S MANUAL. 
EACH R6C0RO CAN THEN BE BR0UGHT INTO MEMORY BY CALLING READ. 

0NE CALL T0 READ BRINGS IN 0N6 RECORD. 

THE EXCEPT 1 0N IN THIS CASE IS DATA READ EXACTLY AS IT APPEARS 0N 
BULK DATA CARDS. F0R THIS* THE PROCEDURE IS CALL PREL0C, THEN 
LOCATE, THEN READ F0R EACH BULK DATA CARD. THIS IS USED F0R IFL. 


DM106100 

DM106200 

DM106300 

DM106400 

DM106530 

DM106600 

DM106700 

DM106800 

DM106900 

DM107000 




DATA 

DATA 

DATA 

DATA 

DATA 


NAM5 i NAM 6 /205, 206/ 
NAM7/NAM8 /207,208/ 

I GR ID / 4501 * 45/ I.FL /10l/ 
TWO* E0R /2* 1/ 

ICON / 102/ IC0NM /l 501 *15/1 BG 


0) 


DM107200 
DM107300 
DM107400 
DM1 07500 

/l 03/ DM 1076 DO 

DM107730 
DM107800 
DM107900 
DM108000 
DM108100 
DM108200 
DM108300 
DM1084DO 
DM108500 
DM108600 
DM108700 
DM108800 
DM 108900 
DM109000 
DM109100 
DM109200 
DM109300 
DM109400 
DM109500 
DM109600 
DM109700 

WR ITE L6 * 601 ) DM109800 

601 FORMAT! 1H I*, - RESTRICTIONS ***-// DM109900 

F- 1. ID N0 . 1 ON GRID CARDS SHOULD BE USED F0R THE CENTRAL RIGID-/ DM 1 1 000 0 
B0DY OTHERWISE THE TRANS FORMAT! 0N MATRICES T AND T- INVERSE -/DM110100 
WILL BE INCORRECT-/ DM110200 

2. THE USE 0F THE T AND T-INVERSE PERMIT THE BASE MOTION WITH -/DM110300 
THE TRANSLATIONAL RIGID BODY M0TI0N SWEPT 0UT FROM THE EQS .-/DM110400 
0F M0TI0N -) DM1 10500 

WRITE 16* 604) DM110600 

604 FORMAT MH * DM110700 

F -3* IT IS ASSUMED THAT THE AXIS 0F ROTATION PASSES THROUGH THE-/ DM110800 

CENTER OF MASS 0F THE ENTIRE VEHICLE. IF THE AXIS 0F R0TAT I 0NDM1 10900 

0 Ml 110.00 

IS T0 PASS THR0UGH A POINT 0THER THAN THE MASS CENTER* THE-/ 
CORRESPONDING GRID P0INT SHOULD BE DEFINED AS SR I D N0. 1-) 

WRITE 16, 605) 

605 FORMAT ( 


THIS IS THE MAIN ROUTINE FOR COMPUTING THE K-PRIME, P-PRIME, G AND 
MASS MATRICES USED IN SOLVING THE ROTATING FLEXIBLE STRUCTURE 
PR 08 LEM 


READ IN OMEGA VALUES AND RG + T-INVERSE FLAGS 


IWW == 104 
I TYPE * i 
jJ = 1 
N = 5 
JINCR - 1 

CALL OPENmOOO* IWW, XKPiLCOL*!), 
FWDREC! $1000* IWW) 

FNAMEt IWW*? R HE AD Cl)) 

UNPACK! $1000* IWW* W* READ) 
CLOSE l I WW* 1 ) 


CALL 

CALL 

CALL 

CALL 


WRITE C0MMENT5 ON OUTPUT LISTING 


F- 

F- 

F- 

F- 

F- 


F- 
F -/ 
F- 
F- 


F- 

F- 

F- 

f- 

F / 

F- 

F- 


4. LOCATION 0F THE CENTER 0F ROTATION (C.R* ) GOVERNS THE-/ 
CALCULATIONS 0F DELTA TRANSPOSE DELTA MATRICES AND THE-/ 
CENTRIFUGAL FORCES WHICH IN TURN AFFECTS THE GEOMETRIC-/ 
STIFFNESS MATRIX! THE OPTION T0 SELECT THE C.R. IS 


DM11 1100 
DM111200 
DM111300 
DM111400 
DM111500 
DM1 11600 
DM1 11700 


AVAILABLE-DMl 11800 
DM1 11900 

THRU DM I CARDS. W(4)*0.0-CG IS NOT CALCULATED AND C.R. IS-/ DM112000 
ASSUMED T 0 BE GRID NO. U W(4) = l.0-CG IS CALCULATED AND C . R. -DM1 1 2100 
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F / 

F- 

F- 

F- 

F- 

F- 


IS C v G • - 

W ( 5 ) *1 * 0 PERMITS THE CALCULATIONS. 0F T AND T-INVERSE.-/ 
tf<5)»0*.0 MAKES T AN.0 T-INVERSE IDENTITY MATRICES*-/ 

IN GENERAL THE FOLLOWING .COMBINATIONS SHOULD BE USED*-/ 
W(4)*W(5>*0.0 0R W ( 45 » W.( 5 ) - 1*0 * W ( 4 ) =W ( 5 ) =0* 0 ASSUMES THE- ) 


DMII2200 
DM112300 
DM112400 
DM112500 
DM112630 
DM112700 
DM112800 
DM112900 

STRUCTURE R0 TATES AB0UT GRID N0. I AND ONLY THE CANTILEVER-/ DM113000 
MODES ARE AVAILABLE. W ( 4) =.W ( 5 ) =1* 0 ASSUMES THE STRUCTURE-/ DM113100 

ROTATES AB0UT THE C.G* WITH TR ANSL AT 1 0NAL RIGID B0DY D0F-/ DM113230 

F- SWEPT 0UT-/ DM113300 

F- 5. IN THE CALCULATION 0F THE CENTER 0.F MASS 0F THE VEHICLE ALL-/DM1 13400 


WRITE 16,899) 
899 F0RMAT ( 

F- 
F- 
F- 


THE GRID P0INTS WITH MASSES ARE USED* IF THERE IS A DUMMY- / DM113500 
GRID P0INT F0R THE DEF IN I T 1 0N 0F THE PLANE F0R CBAR CARDS, -/ DM113600 
D0 N0T PUT ANY MASS AT THAT DUMMY GRID P0INT OTHERWISE WR0NO-DM1 1 3700 

DM113800 

C. MJ WILL BE COMPUTED-) DM113900 

DM l 14000 
DM114100 
DM114200 
DM114330 
DM114400 


F- 
F- 
F- 
F/ 

F- 

WRITE { 6,606) 

606 F0RMAT1 

F- 6* THE NUMBER 0F GRID POINTS IN THE PROBLEM SHOULD BE LESS-/ 

F- THAN 300; THIS ALLOWS UP T0 1800 Q0F *-) 

WRITE (6,602) 

602 FORMAT ( -1H0 , - THE FOLLOWING DATA WAS TAKEN FROM GRID CARDS BY THE DQMU4530 
FUMMY MODULE**:*-/) DM114600 

WR ITE 16 ♦ 612 ) W Cl ) > W ( 2 ) ,W( 3) DM114700 

612 FORMAT! / 1H i -'SPIN RATE VECTOR*** 0ME3A1 = -F10.4, 5X, DM114800 

F -0MEGA2 - -F10.4, 5X, -0MEGA3 = -F10.4) DMU4930 

IF ( W ( 4) *LT* 0.01) WR ITE { 6, 695 ) DM115030 

695 FORMAT 1 1H *-THE FOLLOWING ANALYSIS ASSUMES ROTATION ABOUT GRID N0. DM115100 


i 1.-) 

IF ( W { 4) *GT* 0.0) WRITE (6f696) 


DMI 15230 
DM115300 


C 

C 

C 


FI 

WRITE C6, 611) 

SET UP PACK! COMMON AND LOCATE END OF CORE 


696 FORMAT ( 1H ,.-IN THIS ANALYSIS THE STRUCTURE SPINS ABOUT THE C.G* 3FDM115430 
1 THE SYSTEM DESCRIBED IN THE BULK DATA-) DM115530 

I F t W ( 5 ) * L T . 0.01) WRITE l 6,697) DM115630 

697' FORMAT ( 1H *-THE FOLLOWING ANALYSIS GIVES CANTILEVER MODES 0F A SPIDM115700 
FNNING STRUCTURE BY MAKING T. AND-/1H , -T-INVERSE IOENTITY MAtrUCES-DMl 15800 
F> DM115900 

IF ( W I 5 ) .GT. 0*0) WRI TE ( 6( 698 ) DM116000 

698 F0RMAT ( 1H *-THE FOLLOWING ANALYSIS CALCULATES THE MATRICES T AND TDM116130 
F- INVERSE THUS SWEEPING OUT THE- /IK , -TRANSLATIONAL RIGID B0DY D0F-DMI 16200 

DM116300 
DM116400 
DM116500 
DMi 16600 
DMI 16700 
DM116830 
DM116900 
DM117000 
DM117100 
DM117200 
DM117300 
DM117400 
DM117500 
DM117600 
DM117700 
DM117890 
DM117900 
DMI 18000 
DMI 18100 


611 F0RMAT ( IHO*— N0DE 
132 


D0F 


MASS 


III 


121 


122 


F 131 

TYPIN * i 
TYP0UT = 1 
II * 1 
INCR » 1 

L l = C0RS2 (XXX,* XKP ) 
LC0L » LZ - SYSBUF 
It * LC0L 

IXX = LZ - 2*£YSBUF -2 
I8GR =* IXX - SYSBUF - 1 
D0 750 11*1,3 
750 RGIII) = Q.O 
XM = 0.0 


13 3 


R1 


R2 


R3-/ ) 


229 



on no o noon ooo ooo 


C 

C 

G 


L0CATE FILE THAT GRID P0I.NTS ARE STORED 0N 

CALL PREL0C ! $1000$ XKPML1, IFL ) 

CALL LOCATE! SlOOOfXKPJ HU IGRIO, IFLG) 

L0CATE FILE THAT BASIC GRIO POINT C0-0RDI NATES ARE ST2RE0 0N 


0) 


DM 118200 
DPI 18300 
DM1 18430 
OP 118 500 
DPI 18600 
DPI 18700 
DM118830 
DM118900 
DM119030 
DM119130 
DM119230 
DM119300 
DM119400 
DM119500 
DM119630 
DM119730 
DM119830 
DM119900 
DM120030 
DM120100 

WHEN AN END 0f FILE IS ENCOUNTERED BY FWDREC , READING 0F MG.G ST0PSDM12O2OO 

DM120330 


CALL 0P£NI$1OOO, IBG, XKP C I BGR) « 
call FWDREC! $700, IBG) 

NODE =» 0 

BUILD MASS AND INERTIA TABLE 


UNPACK MGG MATRIX T0 GET MASS AND INERTIAS 
N * 300 

CALL 0 PEN l $1000 , ICONfc XKPMXX), 01 
CALL FWDREC! $1000 V ICON) 

CALL FNAM6UC0N, RHEAD! 1 ) 1 


I C0N) 
IC0NT 


ICON, XMGGf READ) 


950 CALL FWDREC C $699 , 

CALL FWDREC l $699 , 

N0DE 3 N0DE + 1 
NO 1 = (N0DE-1)*6 
CALL UNPACK! $951 
MCN0DE) = XMGGCND1V3) 

XM = XM + M{ NODE ) 

951 CALL UNPACK! $952 V !C0N, XMGG$ READ) 
XMKN0DE) = XMGGCND1+4) 

<952 CALL UNPACK! $953 , ICON, XMGGf READ) 
XC2HN0DE1 •» XMGGIND1+4) 

X I 2 2 ! N 0 D E ) = XMGG1ND1 + 5 ) 

953 CALL UNPACK! $950 fc IC0N, XMGG* READ) 
XI 31! NODE) = XMGG (N01+4) 

XI 32 C N0DE) 3 XMGG!ND1*5) 

XI33CN0DE) = XMGG! NO 1*6) 

G0 T0 950 

699 NODE * 0 

700 C0NTINUE 

CALL READ! $71 0? $710, I.FLf 
NODE * N0DE * I 


A$ 8, 0, IFLG) 


BASIC GRID P0INT INF0* 

CALL READ! $71 Ot* $710, IBGf R, 4, 0, IFLG) 

IF MOTION AT A GRID P0IT IS RESTRAINED IN ALL DIRECTIONS THE L0GIC 
IN THE C0DE CAUSES THE GRID POINT T0 BE DISREGARDED 

WR ITE $6 , 610) ID, IA, MIN0DE), XI1KN0DE), XI 21 ( NODE) , XI22IN0DE) 
W , X I 311 NODE) ♦ XI321N0DE), XI33!N0OE), R(2), R!3), 

W R ( 4) 

610 FORMAT ( II 5 1 118* 10F1GV4) 

IF C W! 4) .GT. 0-0) G0 10 720 
IF I ID .GT* 1) G0 T0 725 
D0 733 11*1,3 
733 RGIII) = R ( 1 1*1 ) 

G0 T0 725 
720 CONTINUE 

D0 701 11*1,3 

701 RGIII) * R G ( II) + M CN0DE )*R 1 1 1 + 1 ) 


DM 1 20400 
DM120500 
DM120630 
DM120700 
DM1 20800 
DM120900 
DM121000 
DM121130 
DM121200 
DM121330 
DM121430 
DM121500 
DM121600 
DM 12 1700 
DM121800 
DM121900 
DM122030 
DM122100 
DM122230 
DM122300 
DM1 22400 
DM122500 
DM122630 
DM122730 
DM122800 
DM122900 
DM i 23000 
DM123100 
DM123200 
DM123330 
DM123430 
DM123530 
DM123600 
DM123700 
DM123800 
DM123900 
DM124030 
DM124130 
DM124230 
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725 C0NTINUE 
G0 T0 700 
710 CONTINUE 
1000 CONTINUE 

I F I W ( 4 ) .EQ. 0.0) G0 T0 719 
00 734 11-1,3 
734 RG(II) * RGC H3/XN 
719 CONTINUE 

WRITE 16, 607) 

607 FORMAT ( 1H „-RG GIVES THE COORDINATES 0F THE CENTER OF R0TATI0N-) 
WRITE 16,615) XM, RGU3,. R GC2H RGI3) 

615 F0RMATC1HO,.- TOTAL MASS = ~,1E15.8, - RGtl) — ,1E15.8, - RGI2) 
F~* 1E15 * 8,< - RGi 3) 1E15.8) 

WRITE 16,620) 

620 FORMAT { 1HG * - M0TI0N CONSTRAINTS IN 1 THRU 6 DIRECT 1 0N AND TOTAL 
F0F AT EACH N0.DE -//■-? N0OE C0NSTRA 1NTS T0TAL DOF -■) 

CALL CL0SEMFL, 1) 

CALL CL0SE ( IC0N, 1) 


0M12433O 
OM 1 24400 
DM 1 24500 
DM124600 
0M1 24700 
DM1 24800 
DM124900 
DM125000 
DM125100 
DM125230 
0ML253D0 
—DM125400 
DM125500 
DM125630 
DDM125700 
DM125830 
DM125930 
DM126000 


**#*♦**» #*#***# ****** ************ if************************************** DM1 26100 


N A Ml 1201 ) 

• it 

NAM8C208) 


REFERS T0 PG-THE FIRST 0UTPUT DATA 8L0CK 3F DUMM0D1. 


REFERS T0 D0F-THE EIGHTH 0UTPUT DATA BL0C< 0F QUMM0D1 . 


THE PR0CEDURE F0R PACKING A MATRIX IS-0PEN THE FILE! EG. 201), 
CALL FNAME, CALL WRITE T0 WRITE THE HEADER RECORD, THEN PACK. 
ALS0 8UIL0 THE TRAILER ARRAY, THEN CALL WRTTRL • SEE 2.2-1 0F THE 
NASTRAN PR 0G R AM MER-S MANUAL F0R TRAILER INF0RMATI0N. 

ALWAYS CL0SE THE FILES 0UT AS S00N AS P0SSI8LE 

CALL CLOSE ( I BGf 13 
CALL Of RECLC0L * N0DE ) 

D0 672 1 1 -1, NODE 

NFREfll) » NDEGF(1, II) + NDEGF C 2, II ) + NDEGF {3,11) 

NFREE - ND8GFI 4, 1 1 ). + NDEGF ( 5(11 ) * NDEGF (6,11) 

672 TDOFUi) - NFREUI) * NFREE 

WR ITEC 6, 625) til, t NDEGF U J , ID , JJ= 1 , 6 ) , TDGF(II), 11*1 , N0DE ) 

625 F0RMAT C 1 15 > 9X, ! 611, 11113 

CALL 0 PE N C$300 ,i NAMi, XKPILC0L + 1H 1) 

CALL FN A MEIN AMI, HEADM) I 

CALL WRI 76 (NAMI , HEAD Cl), TW0', E0R3 

ZER0 0UT R0WS AND C0LS • 0F EACH 6X6 SUBMATRIX WHERE THE DEGREE 0F 


FREED0M IS N0T USED 1 
THE N0N-ZER0 ELEMENTS 


EG F0R NDEGC3) 
ARE THEN M0VED 


= 0 , THIRD R0W 
T0 THE TOP LEFT 


DM126200 
DM126300 
DM126400 
DM126500 
DMl 26630 
DM126700 
DM1268Q0 
DM126900 
DM127030 
DM127100 
DM127200 
DM127300 
DM127400 
DM127500 
DM127600 
DM127700 
0M12780G 
DM127900 
DM128000 
DM128100 
DM128230 
DM128300 
DM128400 
DMl 28500 
AND C0L ARE ZER0DM12863O 
C0 RNER 0F SUB M AT DM 1 2 8 700 



ZER0 0PEN CORE 

DM 128800 
DM128930 


NC0R » 6*N0DE 

DM129000 
DMl 29 LOO 


D0 817 I I-l, NG0R 

DM129200 


XKP(II) * 0.0 

DM129300 

H 

00 

CONTINUE 

. DMl 29400 


D0 818 I I* 1, N0DE 

DM129500 


D0 818 JJ-1,6 

DM129600 

818 

ITT = NDEGF ( JJ,;1 13 + ITT 

DM129700 


CALL KPRIRCN0DE, W) 

DM129800 

300 

CONTINUE 

DM129900 


CALL CL0SEI NAMI » 1) 

DM130000 


CALL 0PENC$5Q1,< NAM 2, XKP$LC0L+1 H 1) 

DM130100 


CALL FNAMECNAM2, HEAD2C1) ) 

DM130200 
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CALL Uft-ITetNAM2t. HEAD2C1H TW0,. E0R) 
CALL 0PEN t $501 / IBG, XKP ( I BGR ).» 0) 
CALL FWDREC ( $501 i IBG) 

CALL PPR I M { RG& N0DE * W) 

501 CONTINUE 

CALL CL0SEINAM2* I) 

CALL CLOSE U 80/ 1) 

ZERO 0UT OPEN CORE 


00 2001 M=a,NI 
01 XKP 111) = 0.0 

CALL 0PENU3OOO* N A M3 i XKP ( LC0L + 1 ) * X) 
CALL FNAME1NAM3* HEA03(1) ) 

CALL WR I TE (NAM3> HEA.03MH T W 0 *« E0R) 
CALL GMATCXM, N0OE, W) 

00 CONTINUE 

CALL CL0SE ( NAM 3 t 1) 

ZERO 0UT 0PEN CORE 


00 3001 M^l, NX 
xKP(i.n = o.o 
CONTINUE 

CALL 0PENt$4OOO, 
CALL FNAME ( NAM4* 
CALL MR I TE { NAM4 * 


HAm\ XKP C LC0L + 1 ) t 1) 
HE AD 4 ( 1 ) ) 

HE AD4 { 1 ) | TW0» E0R) 


CALL AM AT LN0OE ) 

CALL CL 0SE ( N AM4 * 1) 

CALL 0PENt$4OOO,. NAM5 1 XKPCLC0L + 1), 1) 
CALL FNAME ( N AM5 * HE AD 4 ( 1 ) 1 
CALL WRITE ( NAM5 * HE ADM 1 > t TW0„ E0RT 
CALL TMAT(XM,N0DEtW»ITT) 

CALL CLOSE { NAM5 r 1) 

CALL 0PEN t $4000 » NAM6 > XKPCLC0L+1), 1) 
CALL FNAME { NAM6 > HEAD4C l ) 1 
CALL WRITE ( N AM6 *. HEAD4UH TW0» E0R') 
CALL TIMAT(N0DE,WUTT*XM) 

CALL CL0SE { NAM6 ? 1) 

PACK RP MATRIX 


CALL 0PEN 1 $4000 1 NAM74 XKP(LC0L+1), 1) 

CALL FNAME { NAM7 HEAD4(1)1 

CALL WRITE! NAM7 * HEAD4( I ) I TW0„ E0R5 

N1 = NODE 

MMU) * 207 

MM*2) * 0 

MM I 3 ) = N0OE 

MM 14) * 2 

MM (5) * 1 

MM ( 6 ) * N0DE 

D0 3500 1 1*1 » 3 

CALL PACK t RP (1*11)4 NAM7, WRITE* MM) 
3500 C0NTINUE 


PACK MASSES 


OMI 30300 
DM130400 
DM130500 
DM130600 
DM130700 
DM1 30800 
DM130900 
DM131000 
DM1 31 100 
DM1 31200 
DM131300 
DM131400 
DM 1 31500 
DM131600 
DM13 1700 
DM131800 
DM131900 
DM132000 
DM132100 
DM132200 
DM132300 
DM132400 
DM132500 
OM 132600 
DM132700 
DM132800 
DM1 32900 
DM133000 
DM133100 
DM133200 
DM133300 
DM133400 
DM133500 
DM133600 
DM133700 
DM133800 
DM133930 
DM134000 
DM1341Q0 
DM134200 
DM134330 
DM134430 
DM1 34500 
DM134600 
DM134700 
DM134800 
DM134900 
DM135000 
DM135100 
DM135200 
DM1 35300 
DM135400 
DM135500 
DM135600 
DM135700 
DM135800 
DM135900 
DM136000 
DM136100 
DM136200 
DM136300 
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CALL PACKIM, NAM7I WRITE, MM) 

CALL WRTTRL { MMH ) ) 

CALL CL0SE ( NAM7 >. 1) 

PACK O0F 

CALL 0PENt$4OOO*. NAM8I XKP ( LC0L+1 ) , 1) 
CALL FNAME1 NAM8 * HEAD4(1)1 
CALL WRITEINAM8, HEAD4(1){ TW0 , E0R) 

00 3550 II=a,6 
00 3550 Jl=.l,N0OE 
3550 I F ( NDEGF ( 1 1 , J1 ) .GT. 0) O0F(II,J1) = 1.5 
N 1 * 3 
MM ( 1 ) = 208 
MM 12) = 0 
MM (3) * 6 
MM I 4 ) « 2 
MM 15) => 1 
MM ( 6 ) * 6 
00 3600 n=ii,N0oe 

CALL PACK! D0F< U 1 1 HNAM8 „WR I TE , MM) 
3600 CONTINUE 

CALL WRTTRL ( MMl 1 ) ) 

CALL CL0SEINAM8, I) 

4000 CONTINUE 
RETURN 
END 


0M136400 
DMI36500 
DM136630 
DK136700 
DM136800 
DM136900 
OM1370D0 
DM137130 
0M137200 
DM137300 
DM1374C0 
DM137500 
0M137600 
DM137700 
DM137800 
OH137900 
DM138000 
DM138100 
DM1 38200 
DM133300 
DM138400 
DM138500 
DM138600 
DM138700 
DM138800 
OM 133900 
DM139000 


For complete listing of this program on Univac 1108 Computer, 
write to Reference 3, 
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Center of mass in steady state configuration 


Central rigid body 


Center of mass 
at time, t 



Configuration of an appendage 
at time, t 


Steady state configuration 
of the appendage 


FIGURE 1. GEOMETRY OF SPINNING FLEXIBLE APPENDAGE AND CENTRAL BODY 


